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Using Monte Carlo techniques, Ising cubes with ferromag- 
netic nearest-neighbor interactions and enhanced couplings 
between surface spins are studied. In particular, at the sur- 
face transition, the corner magnetization shows non-universal, 
coupling-dependent critical behavior in the thermodynamic 
limit. Results on the critical exponent of the corner magneti- 
zation are compared to previous findings on two-dimensional 
Ising models with three intersecting defect lines. 
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I. INTRODUCTION 

In the thermodynamic limit, critical phenomena may 
occur not only in the bulk of a system, but also at its 
surfaces, edges and corners. To be specific, let us con- 
sider Ising magnets with short-range interactions. Then, 
there are two typical scenarios: (a) bulk, m?,, surface, mi, 
edge, TO2, and corner, TO3, magnetizations may order at 
the same temperature (" ordinary transition" ) , but with 
different power-laws, and (b) surface, edge and corner 
magnetizations may order simultaneously first ("surface 
transition") due to enhanced, strong surface couplings, 
followed by ordering of the bulk magnetization at the 
lower bulk transition temperature ("extraordinary tran- 
sition"). Surface singularities at ordinary and surface 
transitions have been studied extensively, theoretically 
as well as experimentally Q . Most of the rather few 
studies on edge critical behavior dealt with the ordinary 
transition . Only very recently, edge criticality both 
at the surface transition |q] and at the normal transition 
1^ has been investigated. 

Similarly, corner criticality in three-dimensional Ising 
systems has been analysed, to our knowledge, only at 
the ordinary transition, applying mean- field theory 0] 
and Monte Carlo simulations 0. However, that case de- 
serves to be studied at the surface transition as well for 
various reasons. For magnetic properties of nanostruc- 
tured materials, corners are expected to play an impor- 
tant role (l^Jl^ . In addition, magnetic couplings may be 
enhanced at surfaces, especially at step-edges and cor- 
ners Last, but not least, the problem is of genuine 
theoretical interest. At the surface transition, the crit- 
ical fluctuations are essentially two-dimensional. Edges 
are local perturbations acting then like defect lines in 
two-dimensional Ising models leading to interesting non- 
universal critical phenomena |10 ^. Accordingly, corners, 
say, of an Ising cube may be interpreted as intersection 
points of three defect lines. At such points, one also 



expects intriguing non-universal behavior of local quan- 
tities, such as the corner magnetization, following exact 
analytical work on two-dimensional Ising systems with 
intersecting defect lines uM. 



II. MODEL, METHOD AND RESULTS 

We study nearest-neighbor Ising models on simple cu- 
bic lattices with L x M x N spins (usually we shall con- 
sider Ising cubes, i.e. L = M = N) and ferromagnetic 
interactions. The Hamiltonian may be written in the 
form 

— ^ ^ JhSxyzSx'y' z' ^ ^ J s^xyzSx' y' z' 

hulk surface 

^ ^ JesSxyzSx'y' z' ^ ^ JeSxyz^x'y' z' (1) 

edge—surface edge 

where the sums run over bonds between neighboring 
spins, Sxyz = il, with coupling constants to be specified 
below; x{y,z) going from 1 to L{M,N). Free boundary 
conditions hold for the spins in the surface planes. The 
pairs of neighboring spins in the Hamiltonian (1) are lo- 
cated either on edge sites with the edge coupling J^, on 
edge and surface sites coupled by Jes, on surface sites 
with the interaction Jg, or on sites with at least one of 
the spins in the interior of the system interacting with 
the bulk coupling Jf,. We refrained from assigning an- 
other coupling strength to pairs of spins on corner and 
edge sites, which is taken to be equal to Je- 

To study the behavior at the surface transition, T^, we 
chose Js = 2Jb, where ksTs/Jb w 4.975 g, while the 
bulk transition, Tc, occurs at ksTc/Jb ~ 4.5115 Jl5| , p^ . 
The effect of the edge couplings was studied by consid- 
ering the three cases (i) of equal surface couplings, i.e. 
Je ~ Jes — Js, (ii) of reduced edge couplings, especially 
Jes — Js, Je = Jb, and (iii) of reduced edge-surface cou- 
plings, especially Je = Js, Jes = Jb- 

The size of the Ising cubes wih L'^ spins ranged from 
L = 5 to L = 80. In the Monte Carlo simulations, we 
used the efficient single-cluster-flip algorithm. Thermal 
averages were obtained from an ensemble of at least 10^ 
realizations, using different random numbers. In each 
realization, several 10* clusters were taken into account, 
after equilibration. 

The quantity of main interest is the corner magnetiza- 
tion, or more general, the local magnetization mi{x, y, z) 
at site (xyz). mi{x,y,z) may be defined by the correla- 
tion function 
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■mi{x, y, Z) = y/< SxyzSx'y'z' > 



(2) 



where (xyz) and {x'y'z') are topologically equivalent sites 
with maximal separation distance; brackets denote the 
thermal average. In the thermodynamic limit, L — > oo, 
one recovers the standard definition of the local magne- 
tization m(x,y,z) =< Sxyz >■ 'mi{x,y,z) approaches 
closely m{x,y,z) provided the separation distance be- 
tween the two equivalent spins is large compared to the 
correlation length. Certainly, finite-size effects are most 
severe near criticality, as usual. The deviation of mi from 
TO may be monitored by varying the size of the cubes, L, 
and by considering the correlation function between spins 
on equivalent sites with different separation distances (for 
instance, two corner spins may be connected either by an 
edge, by a surface diagonal, or by the bulk diagonal). 

In Fig. 1 the intriguing profile of the local magne- 
tization at the surface, mi{x,y,l), of a 40'^ Ising cube 
with equal surface couplings, case (i), is depicted, at 
ksT/Jt = 4.9, i.e. T w 0.985Ts. The non-monotonic be- 
havior along pathes from the edges or corners towards the 
center of the surface reflects the influence of bulk spins, 
as has been discussed before . Crossover to monotonic- 
ity of the profile, with the largest magnetization at the 
corners, is expected to occur even closer to Tg. At lower 
temperatures, roughly T < Tc, the profile is monotonic as 
well, with the smallest magnetization at the corners due 
to the different coordination numbers at corners, edges, 
and surfaces. 

Near the surface transition, T^, the corner magnetiza- 
tion, say, TO3 = TO/(1, 1, 1), is expected to vanish, in the 
thermodynamic limit, as 7713 cx t^^ , where t is the reduced 
temperature t — \T ~Ts\/Ts. To estimate /Ss, we consider 
the effective exponent [7|P,p^ (3eff{t) — dhim^/ dint. 
When analysing the Monte Carlo data, the derivative is 
replaced by a difference at discrete temperatures. As 
t — > 0, j3eff approaches /Ja, provided finite-size effects 
can be neglected. 

The temperature dependence of the effective exponent 
(ieff for the three different sets of couplings, (i), (ii), 
and (iii), is shown in Fig. 2, displaying only data which 
were checked to be unaffected by finite-size effects. Error 
bars stem from the ensemble averaging performed to de- 
termine m^. The resulting estimates for the asymptotic 
critical exponent (3^ are (i) 0.06±0.01 at Je = Jes = Js = 
2Jb, (ii) 0.14 ± 0.015 at Jes = Js = and Je ^ Jb, 
and (iii) 0.26 ± 0.02 at Je = Js = 2Jh and Jes — Jb- 
In addition, we estimated /Sa at the ordinary transition, 
Je = Jes = Js = Jb to be /Sa = 1.77 ± 0.05, confirming 
and refining our previous estimate based on cornputing 
the corner magnetization from metastable states Er- 
ror bars are inferred from "reasonable" extrapolations of 
the effective exponent, see Fig. 2. 

To explain the Monte Carlo findings on (3^, note that 
the critical fluctuations at the surface transition are es- 
sentially two-dimensional and that corners are intersec- 
tion points of the edges. Now, as had been shown be- 
fore B, at the surface transition edges act as ladder- 



or chain-type defect lines |T^- [20|JTo[ | . The critical expo- 
nent P2 of the edge magnetization is non-universal (being 
non-trivial even in case (i) of equal surface couplings, due 
to the coupling to bulk spins), varying with the edge, Je, 
and edge-surface, Jes, couplings To a given set of 
interactions Je and Jes, one may assign roughly an ef- 
fective defect coupling of ladder- or chain-type, J^'^'^ , 
yielding the same critical exponent for the defect mag- 
netization in the two-dimensional Ising model, /?;, and 
for the edge magnetization at the surface transition of 
the three-dimensional Ising model, (32- Specificly, for a 
ladder-type defect, the critical exponent of the magneti- 
zation in the ladder rows, /?/, is given by 19 20|] 



I3i = 2arctan2(K,"^)/7r2 



(3) 



with Ki tanh(J;/(A:BT2d))/tanh(J/(fcBT2d)), where 
T2d is the transition temperature. Comparing (32 and 
(3i , one may interpret the defect coupling J/ of the two- 
dimensional model as the desired effective coupling J^^^ 
(J is the coupling constant away from the defect line, 
corresponding to Jg in the three-dimensional systems). 

Following this analogy, the critical exponent of the cor- 
ner magnetization, (3^, can be related to that of the mag- 
netization at the intersection of three defect lines in the 
two-dimensional Ising model, (3i, with effective defect 
couplings, J^'^'^. Indeed, in the two-dimensional Ising 
model, the value of Pi has been calculated exactly for 
three intersecting ladder defects by Henkel et al. [Q, 
showing a non-universal behavior, with (3i depending on 
the strength of the defect couphngs, J/. If those cou- 
plings are weaker than in the rest of the system, then 
Pi will increase with decreasing J;(< J), (3i > 1/8, 1/8 
being the well-known Onsager value in the isotropic two- 
dimensional Ising model. In turn, if the defect couplings 
get stronger, then Pi will get smaller. The concrete ex- 
pression for Pi is quite lengthy and will not be repro- 
duced here, but it can be evaluated in a straightforward 
way. 

The effective ladder-type defect couplings J^'^^ in the 
three cases we considered are (i) J^'^^ « 1.22 Jg corre- 
sponding to P2 ~ 0.095 [H in the case of equal sur- 
face couplings, i.e. an effective enhancement of the cou- 
plings at the edges due to the influence of bulk spins, 
(ii) J^^^ « 0.99Js corresponding to P2 w 0.127 § for 
weakened edge couplings, i.e. the enhancement is now 
approximately compensated by the weakening of Je , and 



(iii) J^ 



eff 



0.74Js corresponding to P2 ~ 0.176 M for 



weakened edge-surface couplings, overcompensating the 
enhancement by the reduction in Jes (note that the val- 
ues of P2 differ significantly from those of P3). Using 
these estimates of J^'^'^ , one obtains from the exact ex- 
pression |l^ for the two-dimensional Ising model with 
three intersecting ladder defects of those strengths the 
following values for P, (i) 0.082 , (ii) 0.128, and (iii) 0.21, 
in satisfactory agreement with the Monte Carlo findings 
on P3. Of course, a more refined analysis had to take 
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into account, e.g., the rather comphcated (see also the 
non-monotonic profile in Fig. 1) nature of the edge as 
a simultaneously ladder- and chain-type defect line as 
well as the effect of the bulk spin next to the corner on 
the corner magnetization. Indeed, the good agreement 
between /?3 and (3i in case (ii) may be related to the fact 
that the chain-like character is rather weak in that situa- 
tion. The bulk spin is expected to strengthen the effective 
coupling at the corner especially in case (i), giving rise 
to the reduction in 13^ as compared to (3i. 

Certainly, bulk properties will become critical only at 
the extraordinary transition, at ksTc/Jt ~ 4.5115. For 
instance, the specific heat C is expected to diverge there, 
in the thermodynamic limit. For finite, Ising cubes, 
one observes that a maximum in C near Tc shows up only 
for systems with at least a few thousands spins, getting 
more pronounced as the system size increases (we studied 
case (i) with equal surface couplings and Js = 2Jb). On 
the other hand, the maximum in C near the surface tran- 
sition, Ts, dominates for small cubes, becoming more and 
more suppressed as one increases the size, L. For cubes 
of moderate size, say 15 < L < 60, the temperature de- 
pendence of the specific heat is characterised by an easily 
detectable two-peak structure, with maxima close to Tc 
and Ts.-The height of the two peaks may be easily varied 
by replacing the Ising cubes by slabs. 

In summary, the corner magnetization at the surface 
transition of Ising cubes has been found to display non- 
universal critical behavior, with the critical exponent (3^ 
of the corner magnetization (being distinct from the cor- 
responding edge exponent P2) depending on the strength 
of the edge and edge-surface couplings. The concrete 
value of may be approximated rather well from the 
exactly known value of the critical exponent of the mag- 
netization at the intersection point of three defect lines in 
the two-dimensional Ising model by estimating effective 
defect couplings from the edge critical behavior. 
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FIG. 1. Simulated profile of the local magnetization at the 
surface, mi{x, y, 1), for an Ising model of 40"^ spins with equal 
surface couplings, Je = Jes = Ja = 2Jb, at ksT/Js = 4.9. 
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FIG. 2. Effective exponent /3e// versus reduced tempera- 
ture t for (i) Je = Jes = Js (squares) (ii) Jes = Js, Je = Jb 

(triangles) and (iii) Je = Js, Jes = Jb (circles). Ising cubes 
with up to 80^ spins were simulated, circumventing finite-size 
effects. 
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